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ABSTRACT 


This study investigates the role which mi^t be played by the 
Shuttle in obtaining data which describes the air quality of the 
Eastern United States. A number of aspects of the problem are con- 
sidered including the data requirements of users, a model for 
statistical interpretation of the observations, the influence of 
orbit parameters on the spatial and temporal sampling and an example 
of application of the model. 
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SUMMARY 


The purpose of this report is to present the results of e study 
of the performance of remote sensors of atmospheric pollution when 
placed in Earth orbit. More specifically, the study addressed the 
ability of various instrument /orbit combinations to determine pollu- 
tion levels in the biosphere (herein defined as altitudes from the 
ground to 2 km) of the Northeastern U.S. This particular area is of 
interest because of its extensive air pollution and complex inter- 
urban transport of these pollutants. A considerable national effort 
is underway to achieve a better understanding of pollution transport 
and air chemistry of the region. The addition of spacecraft observa- 
tional capability would add to both studies of historical data (taken 
at the many air monitoring sites within the region) and data obtained 
from experiments specifically designed for regional air pollution 
characterization. Figure 1 illustrates the area of interest. 

The National Aeronautics and Space Administration (NASA), 
sponsor of this study, limited its scope to include specific satel- 
lite orbits and remote sensor payloads. The orbits chosen for study 
were to include those which are consistent with the early missions of 
the Space Shuttle, on which the proposed experiments might be flown. 
The sensors chosen for consideration include two of the more advanced 
atmospheric remote sensors currently under development by NASA. 
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FIGURE 1 

THE AREA OF INTEREST 
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A typical Shuttla misaion appropriata tor thasa axparimants 
might ba of tavan day duration* hava an altituda of 24$ km and an 
inclination with raapect to the plane of tha equator of 57*. The 
reaulting ground track of the orbit ia shown in Figure 2* Based on 
the choice of orbit parameters, tha ground track for thia particular 
orbit will repeat every 111 orbita or, equivalently, after 7 days 
(hereinafter referred to aa the number of days in the repeat cycle of 
the orbit). Figure 2 also illustrates the area covered by the field 
of view as it moves along the orbit track. 

The sensors selected by NASA as examples in this study are 
representative of the state-of-the art in passive remote sensing. 

The selection of specific instruments was related to their state of 
development and ability to detect carbon monoxide (CO) and other 
important air pollutants. Carbon monoxide was chosen as the most 
important indicator of sensor capability because of its importance in 
the air quality of the Northeastern U.S. 

After a review of the possible candidate sensors, two were 
chosen for review in this study. One of the sensors, MAPS (Measure- 
ment of Air Pollution from Satellites), is designed to utilise the 
infrared radiation emitted by the Earth and its atmosphere to detect 
trace quantities of carbon monoxide. The principle of MAPS operation 
is the passage of the atmospheric infrared radiation spectrum through 
a sample of carbon monoxide in the instrument. The sample selective- 
ly absorbs that part of the energy which was emitted by carbon 
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■onoxid* in tht atBotpharni Chtrnby providing a ataiurt of tht total 
aaount of that gat within tha viaw of tha inatruMnt. 



Tha other tentor, CIMATS (Correlation Inter ferooM ter for the 
Heaeureaent of Ataoephe^ric Trace Speciea)^ obtains inforaation on the 
radiation upwelling froa the Earth and its ataosphere in the fora of 
an interferograa (Fourier transfora of the radiation speetrtn). The 
interferograa, which covers only a selected portion of the infrared 
epectrua, is then coapared with intarferograaui recorded under labora- 
tory conditions in order to deteraine pollution levels. CIMATS is 
designed to sieasure the trace gases CO, aethane (CH4), nitrous oxide 
(N2 0), aaaonia (NH4) and sulfur dioxide (8O2). 

The MAPS instrusMnt has a noainal full angle field-of-view of 
4.33*. CIMATS has a noainal field-of-view of 7*. In both cases the 
field-of-view can be reducedto about 2* with the aid of e telescope. 
Froa a noainal orbit altitude of 400 ka, these f ields-of-view trans- 
late into a 14 ka (for 2*) to 50 kn (for 7*) resolution eleaent on the 
surface of the Earth. 

The sensors chosen for consideration are soaewhat liaited in 
capability as a result of their basic operating principles. A nwaber 
of these liaitations apply to both instruaenta regardless of the fact 
that the sensing techniques rely on considerably different concepts 
and/or parts of the electroaagnetic epectrua. 

Neither instruaent is capable of directly observing the vertical 
variation in pollution Uv*ls. The detected radiation, which is 
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indicative of the level of {Mllution, eaenetes with an efficiency 
which varies with altitude, each voltme element of the atmosphere 
contributing to some degree. Each measureMnt, then, is the integral 
over altitude of the pollution level, weighted by the sensitivity of 
the instrument to each altitude. The pollution concentration inte** 
grated over altitude is called the column burden. 

The amount of infrared radiation reaching a sensor in orbit is 
determined by the sisniltaneous absorption and emission of radiation 
by the constituents of the atmosphere. Riotons emitted near the 
ground have a highly variable probability of being absorbed before 
reaching the sensor, depend* ng on the wavelength and the concentra- 
tion of atmospheric constituents which absorb at those wavelengths. 
Photons emitted higher in the atmosphere have a higher probability of 
reaching the sensor both because Che path is shorter and because the 
density of Che atmosphere (and its pollutants) decreases rapidly with 
altitude. When these effects are translated into graphic form 
(called a weighting function) as illustrated in Figure 3, it can 
be seen that infrared instruments are generally insensitive to 
pollution near the ground and may exhibit maximum sensitivity at 
altitudes in the 0.5 to 2.5 km range. 

The curve illustrates that the altitudes of interest (0-2 km), 
as well as many others, contribute to the signal detected at Che 
spacecraft. More important, Che contribution from the very low alti- 
tudes, which are often most polluted, is negligibly small compared to 
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FIGURE 3 

WEIGHTING FUNCTION FOR 0.1 ppm CO UNIFORMLY DISTRIBUTED 
IN THE TROPOSPHERE FOR A LAPSE RATE OF 2 K°/km AND A 
SURFACE TEMPERATURE OF 298 K. 
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that of higher altitudes. While the curve is only an example for a 
particular pollutant at a particular wavelength, it is typical of the 
sensitivity of infrared instruments. 

As an alternative, the CIMATS instrument can utilize reflected 
solar radiation rather than infrared radiation. In order to have 
adequate radiation intensity to provide a useful measurement, the 
solar elevation relative to the point directly below the sensor must 
exceed approximately 30 degrees. This factor limits the usefulness 
of CIMATS. Both orbit characteristics and the seasonal variation in 
solar zenith angle determine the local solar elevation. Figure 4 
shows the seasonal variation in elevation of the sun at its zenith 
for a typical orbit. The figure illustrates t’*at observations are 
difficult to obtain unless the date of the mission is chosen with 
care, and also that the length of a mission utilizing CIMATS is 
limited by solar elevation effects to about 20 days. Measurements 
derived from CIMATS using solar radiation are also limited in their 
ability to provide a measure of biospheric pollution. The column 
burden detected by CIMATS includes a contribution from all altitudes, 
the weighting function depending on the atmospheric pressure and the 
relative pollution concentration. 

The inability of the instruments to make direct observations of 
particular altitudes or ranges of altitudes means that pollution levels 
near the ground, which are of particular interest in this study, can 
only be obtained after processing of the data. Processing methods 
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FIGURE 4 

TITUDES COVERED BY A TYPICAL 56° ORBIT WITH THE REQUIREMENT 
THAT THE SOLAR ELEVATION EXCEEDS 30 




which mi^ht provide such estimates of pollution levels have not yet 
been developed, although it is theoretically possible to obtain such 
results* 

Roth instruments would require auxiliary data in order to obta^^ 
even column burden intormation* For example, information on the 
presence ot clouds in the tie Id of view, their location, and altitude 
are essential since the range of the atmosphere over which the instru* 
ments obtain information is interrupted by clouds. Other information 
necessary for data processing includes temperature and water vapor 
profiles of the atmosphere. 

Yet another limitation of both instruments is that neither 
lias been flown in space. There is doubt about the quality ot the 
column burden data which might be i'^btained. This lack ot t light 
experience also means that there is no data base which can be used to 
establish data processing methods which result in the best possible 
estimates of pollution levels. 

In order to quantify the performance which might bt' expt'cted 
from various sensor/orbii c omh i nat i ons , a system model was dt'veloped 
during this studv. It mathematically describes the potential per- 
tormance bv including sampling cha rac t e r i s t i c s (as dt' fined by o,bit 
paranu't^'rs and thr inciilence of c 1 t>ud i nes s ) , instrument tactors (such 
as t i t' I d -o f -V i ew and t inu' constants') and the statistics ot pollution 
variation which ari' typical tor the region of intt'rest. The result- 
ing application ot the modt' 1 is a calculatum the uncertainty 
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associated with pol hit ion estimates made under various observa t ional 
conditions. In this context uncertainty refers to an estimate of the 
measurement error. Direct measurement of the error is not possible 
in this experiment, since the true value cannot be known, and an a 
priori estimate of the error must suffice. In utilizing the model, it 
was assumed that appropriate data reduction techniques cald be 
developed so that pollution levels could be detected for the lower 
altitudes of tlie atmosphere. Obviously, in the case that such data 
analysis techniques are not developed, such missions would be of 
considerably lt‘ss value. As is commonly the case in remote sensing 
programs, an effective data analysis procedure is essential. 

The sensor/orbit model takes into account the sampling which 
can be achieved (as determined by orbit conditions and the influence 
of cloudiness). The sampling is charac t e ri zed by the spatial distri- 
but ion of observations over the area of interest as well as the local 
time of the samples. Temporal sampling is delined by the frequency 
with which various areas are observed, the local sun time of the 
-''hserva t ions , and the solar elevalivMi during the obse rva t i ons . 
vSampling for various t'rbils was derived trom an orbit propagation 
c a 1 cu I a t ion. 

, Kt'sults of i lu' nuHlel appt'ar in Table I where t lU a I uncertainties 

(resulting t lemi all sources) art' d ocunu'tU t*d . Tlu' Table illustratt's 

I the total unct'rtaintv in pollution estimatt's as a t unction oi mission 

durat iv>n and re[U‘at cvclt' for a st't t'tt pi'*! hit ion statistics dt'rived 

I 

I 
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TABLE 4-111 


TOTAL UNCERTAINTY (IN PERCENT) IN POLLUTION ESTIMATES 
ASSUMINC 50% CHANCE OF CLOUD COVER AND 50% UNCERTAINTY 
IN INDIVIDUAL POLLUTION ESTIMATES 


MISSION DURATION (clays) 


REPEAT CYCLE (days) 

7 

15 

30 

L 

27.9 

19.4 

13.8 


J7. I 

25.5 

18. 1 

3 

45.1 

JO, 9 

21.8 

•4 

52. b 

,35. b 

25.2 

s 

56.0 

39.8 

28. 1 

b 

59.2 

4 h 5 

30.8 

7 

62. J 

47.0 

33. 3 
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from real data* The analysis includes the effects of incomplete 
sampling, instrument noise and delays in instrument response to 
spatial and temporal variations in pollution* The results are 
presented assuming a 30% uncertainty in individual measurements* 

This value may vary from instrument to instrument, and the values of 
Table I may change accordingly. However, the Table clearly illu- 
strates the expected trends, particularly that missions longer than 
those currently planned are required to obtain the lowest possible 
uncertainty in pollution measurements* Another important factor is 
the trade-off which exists between the density of coverage of the 
area of interest and the data quality* The shorter repeat cycles 
(which result in Lne lowest uncertainties) also exhibit the poorest 
density of sampling of the region of interest. Choosing an orbit 
which is characterized by a high density of spatial sampling will 
result in relatively high uncertainty in pollution estimates* 

The value of the sensors and missions considered in thir study 
can only be qualitatively judged after comparison with performance 
needs associated with the regional air pollution problem. Three 
categories of data uses can be defined, each with its own set of 
requirements* The first, and most demanding, of these is scientific 
applications* The data must be of high quality (low uncertainty), 
exi)ibit spatial and temporal sampling which allows a complete descrip- 
tion of pollution variability in space and time, and have sufficient 
duration to satisfy the needs of the particular experiment. The 
second general area of application of the data is monitoring. The 
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primary needs are that the observational period be very long, the 
measurements of sufficient quality to determine if standards have been 
exceeded and the spatial sampling pattern representative (though not 
necessarily highly dense). Third, public information data needs 
require only that a general description of pollution location and 
movement can be made. Missions only need be long enough to describe 
the movement of pollutants over inter-urban distances. Quantitative 
concentration data is not essential but the data must supply a general 
characterization of the problem. 

In view of the performance which can be expected from the sensor/ 
orbit combinations under consideration, only the public information 
role can be satisfactorily performed. Scientific applications 
require data quality which cannot be achieved by the instrument/orbit 
combinations, particularly when compared with the capabilities of 
ground-based sensing systems. Monitoring needs cannot be met because 
of the rather short duration of the spacecraft missions. 

The ins trument/orbit combinations do, however, hold the possi- 
bility for providing qualitative data on the movement of polluted air 
masses. Under the conditions imposed by the sponsor on instrument 
selection and orbit characteristics, this qualitative data will be 
sufficient to inform the public of the long-range transport of 
polluted air, illustrate the general relationship between meteorology 
and air pollution and, possibly, identify major pollution sources. 

The results of the study also show that mission durations which 
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greatly exceed the current proposed seven day flight, will result in 
data quality which, while not being sufficient for scientific or 
regulatory purposes, may be adequate to provide a reliable picture of 
pollut ion levels. 

The value of the data to the public will, of course, be a 
function of its accuracy and the fidelity with which it describes the 
region of interest. It is evident, then, that data analysis proce- 
dures which allow determination of pollution levels near the ground 
are greatly needed. The experiment as described in this report is 
arbitarily limited by the fact that only nadir--viewing instruments are 
considered. The inclusion of some method of cross-track scanning 
would allow the selection of short repeat cycle orbits (which minimize 
uncertainty in the pollution estimates) but still provide good spatial 
coverage of the region of interest by scanning the f ield-of-view to 
the areas which the satellite does not fly directly over. This 
approach is not a complete solution to the spatial-temporal sampling 
problem, however, since the number of samples is still limited by the 
duration of the mission. Scanning off the orbit track will reduce 
the number of observations of areas directly below the spacecraft in 
order to provide more complete spatial sampling between the orbit 
tracks. The most appropriate solution to limited sampling is increas- 
ing the duration of the mission. 
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l.O INTRODUCTION 



Providing a clean environment requires that consideration be 
given to the sources, transport, transformation and sinks of the 
various pollutants. This is particularly true the case of the 
atmosphere, which is a repository for many of man's pollutants 
(called primary pollutants) as well as being a medium in which 
chemical and physical changes can occur, resulting in secondary 
pollutants. The atmosphere has become a sort of chemical factory. 
During a stay in the atmosphere, various pollutants may undergo 
photochemical changes, evolve into particles, dissolve in water 
droplets, become oxidized or reduced, etc. Ultimately, they may find 
their way into populated areas, create "smog," and cause health 
effects, material degradation, etc. The complexities and unpredict- 
ability of these processes are increased by the fact that the source 
and the affected area may be far apart due to transport and trans- 
formation. The effects are seen in at least two ways; 1) the totality 
of all sources ccxnbiiie to create a minimum pollution level or "back- 
ground" which is uniform over large areas (globally in some cases but 
more commonly over a regional area) and varies slowly with time, 2) a 
specific source or group of sources produce a plume of pollution 
whose path can be followed over long distances and which may eventual- 
ly pollute another area. The distances involved may include the 
impact of a power plant on areas a few kilometers away to inter-urban 
transport over a range of hundreds of kilometers. The inter-urban 
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transport of pollutants is referrc.i to as rcfcional pollution and is 
the subject of this report. Considerable attention has been given 
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recently to the problesi of regional pollution. A survey is included 
in a recent MITRE Report 11). The report notes that the conventional 
obswivation, modeling and control strategies are less reliable on a 
regional (inter-urban) scale than when applied to a single urban 
area, where the scale of the problem is not too large to be dealt 
with effectively. 

Air quality modeling and prediction becomes more difficult since 
the norsul assumptions used in urban air quality models no longer 
apply: 

e Complex and generally unknown conditions exist for thu 

important meteorological coruitions of wind, precipitation 
and atmospheric stability. 

e The number of possible fates of the pollutants must be 
expanded to include deposition on various surfaces (lakes, 
rivers, forests, etc.), washout in rainfall or snowfall and 
chemical or physical transformations which normally cannot 
occur in the time scale of pollution transport across an 
urban area. 

• Typical diurnal variations in wind, temperature and atmo- 
spheric turbulence are different for urban and rural areas. 
Daily (diurnal) variations in wind or atmospheric turbu- 
lence must be considered. Other atmospheric conditions 
also exhibit variations between urban and rural areas making 
their influence on pollution chemistry difficult to assess. 

e Natural sources of air pollution may be of significance 
at the scale of space and time characteristic of regional 
air pollution. 

Control of air pollution produced in another city, state or 
region of the United States clearly invites the participation of 
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regional or Federal authorities. They are faced with the difficulty, 
however, that the complex transformations which may occur over some 
indeterminate route produce secondary pollutants whose source cannot 
be identified. This problem can be resolved to a great extent, of 
course, by observational programs which follow the "plume** of an 
urban area or industrial complex, monitoring the chemical and phys- 
ical transformations which may occur in its many pollutants, and 
determine those localities where the pollution has an impact on air 
quality. Even pollutants emitted from tall stacks, as is common for 
power plants, may have a distant effect through mixing in the atmo- 
sphere, settling (in the case of particulates) or changes in vertical 
components of the wind. A number of studies of these problems have 
been carried out and are cited in reference (1). 

There are clear limitations, of course, in the ability of an 
observational effort to produce an adequate characterization of the 
processes which result in pollution remote from the source. A 
proper experiment would have to include: continuous monitoring over 

the entire area of interest, monitoring of vertical and horizontal 
variations in space and time of all pollutants and their possible 
derivatives, and determination of factors which control chemical 
conversion including catalytic materials, meteorological conditions, 
solar intensity, etc. Even the most ambitious effort cannot collect 
all of this data on a continuous basis, particularly if conventional 
grcund-based sensing systems are used alone. For this reason, most 
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recent regional air pollution ttudiea have included at least some 
airborne sensing for improved spatial coverage although these air* 
borne sensing systems conmionly use the same instruments as used in 
ground*based observations. These instruments are either ^'contact 
sensors'*, which rely on moving a gas sample through a chemical or 
optical sensing system or "grab samples" which require that a gas 
sample be obtained, stored and brought to a laboratory device for 
analysis. Experiments of this type are limited in their duration 
either because of periods of bad weather, instrument maintenance, or 
more commonly, lack of funds for aircraft, personnel, and processing 
of large quantities of data. They do provide advantages over ground* 
based sensor systems since three dimensions of pollution variation 
can be detected and large areas can be quickly covered with ^ach 
set of sensors. 

Another approach to regional air pollution studies is the use 
of historical data assembled f^om selected ground-based observation 
stations distributed over the entire region. These stations are 
usually in place as part of the air monitoring program of a partic- 
ular city or state and have been independently calibrated so that 
their data may not be entirely consistent. They do provide hoaever, 
a record of pollution levels uver long periods of time, something 
not provided by most regional experiments. There may also be con- 
siderable limitation in the type of primary and secondary pollutants 
which are sampled. For example, a coimsunity may want to monitor sul- 
fur dioxide (SO^) to determine if it meets state or Federal standards, 
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with little concern for the sulfate aerosols which evolve from SO^ 
(since that transformation may occur many kil^imeters away, many hours 
later)* However, that same community may be polluted by sulfates 
which began as SO^ far away at an earlier time* A complete experi- 
mental program would attempt to monitor both the upwind and downwind 

SC^ and sulfates and the rate of transformation of SO^ to sulfate in 
2 2 

space and time* Clearly each approach to regional air pollution 
studies is limited and invites improvement, basically in the area of 
spatial and temporal sampling* 

This study is designed to provide consideration of yet another 
means of contributing to both historical studies and specific experi- 
ments by making use of remote sensing techniques* Airborne remote 
sensing has been utilized in some regional experiments but is rela- 
tively uncommon (contact sampling or grab sampling for laboratory 
analysis is much more common). Satellite remote sensing has sup- 
ported such experiments by providing information on clouds, wind 
fields, snow and ice cover, etc. The remote sensing techniques 
discussed in this report can potentially provide information on 
pollution concentrations by making use of the radiation upwelling 
from the atmosphere and the surface of the Earth (both infrared and 
visible wavelengths can be used)* The detection is performed by 
spectroscopic analysis of the radiation and can provide data on the 
fraction of the atmosphere which is ma-^e up of specific pollutant 
gases. The principles of operation for each instrument are discussed 
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in Section 2, along with data on the gases which can be detected, 
sensitivity, etc. The spatial and temporal sampling which can be 
achieved in applying these instruments will be a function of the 
satellite orbits used. Section 3 provides a discussion of the 
spatial and temporal sampling which can be expected, with specific 
reference to an area of the U.S. ^ich is commonly affected by 
inter-urban transport of air pollutants, the Northeast. 

Characterization of the data quality which can be expected from 
the instrument /orbit combinations results from a statistical data 
interpretation model, described in Section 4. It describes the 
uncertainty associated with estimates of the pollution as a result 
of three contributing factors: 

• uncertainties associated with incomplete spatial and 
temporal sampling 

• uncertainties in individual pollution estimates 
associated with random and bias inscrument errors 

• uncertainties associated with the response of the 
instruments to the spatial and temporal variability 
in the natural pollution concentration 

Section 5 discusses the data needs associated with various types 
of regional air pollution studies. Three possible applications are 
considered: 1) scientific studies intended to improve our under- 

standing of the physical and chemical processes involved in regional 
air pollution, 2) monitoring studies intended to better characterize 
the sources, transformation and fate of air pollutants and 3) public 
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information of the data as a tool in raising the public's awareness 
of the long range impact of air polluti^nts* 

Section 6 presents a comparison of the data needs of Section 5 
and the expected data quality as derived in Section 4. 
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2 . 0 INSTRUMENT PERFORMANCE 


The acquisition of meaningful air quality data will require that 
the instruments chosen for the missions can detect the appropriate 
gases, have relatively low noise and sensitivity to interfering 
species, and can differentiate the low altitude component of the 
column burden. It is the intent of this section to characterize the 
anticipated operating capability of the sensors under consideration 
(CIMATS and MAPS)^ as well as discuss their ability to provide data 
on the pollution of the biosphere. This information will be used 
in Section 4 as input to the system performance model. 

2. 1 General 

Over the past several years many instruments have been developed 
for remote sensing from satellite platforms. Previous reports 
have treated these sensors and their on^rating characteristics in 
great detail [1,2]. Given the interest in carbon monoxide as an 
urban and regional pollutant and the relative state of development of 
the various sensors suitable for space missions, the choice of 
instruments to consider for the purposes of this report was con- 
siderably narrowed. Only MAPS and CIMATS will be treated in this 
study as they are, at present, the only passive instruments with a 
designed capability for the measurement of CO. Each has also had 
sufficient field testing to allow some degree of confidence about the 
performance which might be expected. 

♦Correlation I^nterferometer for the Measurement of Atmospheric Trace 
Species; Measurement of Air Pollution from Satellites 


2.2 MAPS 


The MAPS instrument uses gas filter correlation in order to 
obtain high radiation throughput and resolution. A sample of the gas 
of interest » say CO, is placed in a gas cell located in an optical 
path carrying radiation from the scene to the detector. Another, 
otherwise identical, path traverses an evacuated cell. The difference 
in the electrical outputs of the detectors is the signal difference 
which provides an indication of the amount of the gas of interest 
between the infrared radiation source (normally the Earth and its 
atmosphere) and the instrument. A further refinement in the MAPS 
instrument involves the inclusion of a gas cell at a different 
partial pressure of the gas of interest than is contained in the 
first gas cell. By comparing the difference signals derived from 
each gas cell and the evacuated cell, the measurement is less subject 
to errors caused by the presence of interfering species in the 
atmosphere [3]* An added potential advantage in the use of two gas 
cells of different pressures is the ability to obtain some coarse 
information on several vertical layers of the atmospheric path. 

Figure 2*1 shows the optical layout of the MAPS instrument. 

The reference blackbody is sampled alternately with the scene radia* 
tion to provide a baseline reference level. The hot and cold black- 
bodies are adjusted to span the anticipated temperature range of the 
scene* They are sampled at a higher frequency than the scene energy 
and reference blackbody in order to allow electronic separation of 
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FIGURE 2-1 

MAPS OPTICAL LAYOOT151 


the signals for subsequent processing. Gain controls on the elec- 
tronics are driven by the balance sources so as to keep the gain of 
the gas cell paths equal to that of the vacuum path for both temper- 
ature extremes. This procedure effectively balances the system and 
reduces the effects of variations of radiance reaching the two 
detectors I4l. For the optical path through the evacuated cell, the 
modulated energy from the balance sources is used for radiometric 
calibration, while that from the scene and reference is used to 
determine scene brightness. 

The MAPS instrument is scheduled to fly on a Shuttle mission in 
1979. Table 2-1 shows the instrument characteristics which will be 
used for that experiment. In recent aircraft tests, the two gas 
cells have been filled with CO at partial pressures of 0.35 atmos- 
pheres and 0.1 atmospheres. Current efforts are directed at deter- 
mining the optimum pressures for these two cells. 

In the nadir-viewing (or vertically down- looking) mode, MAPS 
will provide data on the total column burden of CO in cloud-free 
areas of the Earth. Planned refinements in the choice of gas partial 
pressures and data processing may result in data representative of 
two broad vertical bands in the atmosphere. The importance of these 
developments is discussed in Section 2.4, 

2 . 3 Cl MATS 

The CIMATS instrument is a two channel infrared (IR) interferometer 
which uses the same general principles of operation as its predecessor, 
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TABLE 2-1 


MAPS EXPERIMENT CHARACTERISTICS [e] 


Fi eld-o f “ vi ew , de gree s 

Spectral range , micrometers (pm) 

Aperture, centimeters 

Sensitivity, (NEN) , watts/ 
steradian-centimeter2 

Detectors (3) 

Detector temperature, degrees 
kelvin 

Detector cooling 
Response time, seconds 
Data sampling interval, seconds 
Signal channels 
Temperatures 

Data storage, bits per second 
(self-contained tape recorder) 

Input power, watts 

Standby 

Balance/Calibrate 

Operate 

Size, Including single baseplate, 
cm 

Weight, kilograms 


4.33 
4. 5-4. 8 
2.94 X 2.94 

1.7 X 10~® 

PbSe 

195 

Thermoelectric 
8 to 10 

1.0 

15 

50 

2814 

80 

95 

105 

75 X 75 X 45 
65 
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the COPE* instrument (7,8,9j. The A channel uses reflected solar 

energy in the near-IR (2-2.4 fjim) and the B channel uses thermal 

radiation of the Earth and its atmosphere (4-9 m). Table 2-II 

indicates the various species for which CINATS is currently being 

built and the filter parameters associated with each measurement. 

The A channel will use lead sulfide (PbS) detectors cooled to 195®K 

-13 

with a noise equivalent power (NEP) of 5 x 10 watts per root hertz 
(w/vHz). The B channel uses mercury - cadmium - telluride (HgCdTe) 
detectors, cooled to 77®K by liquid nitrogen, having a NEP of 3.5 x 
10 w/yAh 2 at 8.9 m. The instrument f ield-of-view is nominally 7® 
with a 2® f ield-of-*view option when a second telescope objective is 
used • 

The instrument is a variant of a normal Michelson interferometer. 
A time-varying path difference is introduced into one of the optical 
paths by the rotation of a phase plate. The other optical path is of 
constant length (phase). Recombination of the two paths results in 
an alternating reinforcement and cancellation, depending upon the 
phase difference cf the two rays. The result is analogous to a 
Fourier transform of the incident radiation. The effect of this 
design is to transform the electromagnetic variation at optical fre- 
quencies to electrical signals which can be conveniently processed. 

The basic components of the instrument are illustrated in Figure 2-2. 
The range of scan (phase delay) is chosen to optimize the signal-to- 
noise ratio for a particular species of pollutant in the atmosphere. 
*Carbon Monoxide (CO) Pollution Experiment. 
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In the presence of interfering species , the interferogrem 
does not provide a unique indication of the species of interest* 

Table 2-1II lists the major inter ferent and the expected signal-to- 
noise ratio for each measurable species* To reduce interference 
effects, use must be made of weighting functions which are chosen to 
minimize the contribution made ^ / constituents other than those of 
interest, while maximizing the signal-to**noise ratio of the desired 
species* The weighting functions are determined both theoretically 
and experimentally* The best available atmospheric models and 
spectral absorption properties are used to construct a "base inter* 
ferogram." Variation of the concentrations of various constituents 
results in a set of different interferograms* Since each inter- 
ferogram consists of a set of data points, the set of interferograms 
may be considered a set of simultaneous linear equations which, in 
turn, may be solved for the appropriate weighting functions* Labora- 
tory experiments are then performed to generate interferograms based 
upon actual gas concentrations at a range of partial pressures* 

In operation, interference filters limit the range of the 
radiation incoming to the instrument to a very narrow band, selected 
to coincide with the region where the effects of absorption by the 
species of interest are greatest* Interferograms are generated while 
observing the area of interest in a nadir-viewing (vertical down-look- 
ing) mode, one complete scan requiring 1 to 3 seconds* These inter- 
ferograms are compared to the stored weighting functions to determine 
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TABLE 2-111 


INTERFERANTS ANTICIPATED SIOJAL-TO-WISE 
RATIOS [s/n] for CIMATS MEASUREMENTS 110 } 


Channel 

Species 

Assumed 

Background 

Center 

Wavelength 

Interferants 

S/N 

A 

CO 

0.1 ppm 

2.33 vm 

CH^, H2O 

20:1 

A 

CH^ 

1.5 ppm 

2.33 ym 

CO2, H2O 

40:1 

A 

N2O 

0.2 ppm 

2.27 um 

CH^, NH3 

10:1 





CO2, H2O 


A 

NHj 

1.0 ppb 

2.23 ym 

CH^, H2O 

2:1 

B* 

CO 

0.1 ppm 

4.62 yo 

CO2, H2O 

2:1 

B* 

N2O 

0.2 ppm 

4.62 ym 

CO2, H2O 

2:1 


SO2 

2.0 ppb 

8. 70 ym 

CH^. N2O 

not 

aval 1- 





O3, H2O 

able 


*S/N ratios of thermal channel are based on the assumption 
Independent measurement of temperature profile of the atmosphere. 
Availability of such measurements will increase S/N ratios. 
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the best fit of the data. Scale factors may be applied to the 
weighting functions in order to facilitate this matching process. 

When the correct weighting function is applied to the data, the 
result gives the column burden of the specie of interest. In some 
cases, CIMATS will require vertical profiles of temperature, pressure 
and humidity to perform calibrations. 

2.4 Current Limitations of the Instruments 

Both MAPS and COPE have been used for the measurement of the 
column burden of CO from aircraft platforms but neither has yet been 
flown on spacecraft. MAPS is tentatively scheduled for inclusion on 
Orbital Flight Test (OFT) 2 of the Space Shuttle in July 1979. No 
data base, therefore, exists which could serve to indicate the 
accuracies of these instruments when used on a satellite. Based upon 
the experience of each instrument at aircraft altitudes, there is 
little doubt that each can obtain measurements which are related to 
the total column burden of the gas of interest. 

For the purposes of measuring air pollution as it affects the 
biosphere, however, column burdens may not be sufficient. As a 
minimum, the concentrations of various pollutants within the biosphere 
will be required, and for the purposes of public health, concentra- 
tions at lower altitudes will be of g7;eater importance. Tlie ability 
of the instruments to measure pollution levels at low altitudes 
depends on their operating characteristics. For example, the use of 
infrared radiation emitted by the Earth and its atmosphere results in 
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a vertical sensitivity or "lighting function" as shorn in Figure 
2-3. The maxinum sensitivity occurs for altitudes ranging from 0*3 to 
2*3 km, but the instrument is insensitive to pollution very near the 
ground because the gas temperatures are nearly equal to the ground 
temperature. The shape of the curve results from the fact that the 
amount of infrared radiation reaching the sensor in orbit is deter- 
mined by simultaneous absorption and re-emission by the constituents 
of the atmosphere, particularly the pollutant gas. Radiation emitted 
near the surface of •‘he Earth has a relatively low probability of 
reaching the sensor without first being absorbed, since the radiation 
for which the sensor is tuned is, by definition, that part spectrum 
which is characteristic of the gas of interest. The radiation emitted 
from higher altitudes has increased probability of reaching the 
detector, both because it traverses a shorter distance and because the 
influence of the pollutant gas may be reduced for two reasons: 1) its 

fraction of the atmosphere may be declining because it is concentrated 
near the ground, and 2) the density of the atmosphere as a whole 
decreases with altitude. Tlie result of these effects results in a 
curve typified by Figure 2-3 (if the assumptions indicated in the 
figure caption apply). The signal produced by the sensor will be a 
sum of th*' signal levels contributed by each layer of the atmosphere. 
Because the actual conditions of the atmosphere during the observa- 
tions are not known, the sensor signal cannot directly give informa- 
tion on the pollution level of the near-ground region. Use of solar 
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ALTITUDE (km) 




radiation, as can be done with CIMATS, also results in a signal which 
is not totally representative of the pollution level near the ground* 
The radiation reaching the sensor has passed through the atmosphere to 
the ground, back through the atmosphere to the sensor being absorbed 
by the pollutant(s) of interest along the way* Because the rate of 
absorption per meter along the radiation path is determined by the 
number of pollutant molecules, absorption of radiation is a function 
of both the atmospheric pressure and the relative pollution concentra- 
tion* Again, these conditions are not known in general so that 
special data processing techniques will have to be developed to 
provide near-ground estimates* 

While no claims are presently made that either instrument can 
obtain vertical profiles of pollutants, it may be useful to consider 
the possible methods which could be used to extract some vertical 
information frcm either MAPS or CIMATS. 

For both MAPS and CIMATS, clouds pose a problem in coverage* 

For reflected solar radiation, they, rather than the Earth, provide 
the reflecting target. For emitted thermal radiation, clouds appear 
opaque and form the lowest altitude from which radiation reaches the 
instrument. It might be possible to convert this apparent liability 
into an advantage. It is possible to obtain some information on the 
concentration of pollutants in the layer between the Earth*s surface 
and the cloud cover by taking the difference in signals between 
scenes of complete cloudiness and complete lack of clouds. 


-ffrii-iff^riiitfir ‘ Ilf 




In order for this system to be useful, some auxilary information 
would be required* One would have to know cloud locations within the 
f ield-of-*view. Cloud heights, temperatures and extents would also be 
required to effectively analyze the data. This technique is, of 
course, limited by the meteorology which happens to occur at the time 
of the overflight and is, therefore, unpredictable. The use of a 
scanning small f ield-of-view instrument would help in this regard 
since the instrument could be oriented to the clear or cloudy regions 
in a scene. 

For MAPS, another possibility for vertical information of 
gas concentration may be found in the partial pressures with which 
the gas cells are filled. Since the gas cells provide the filter 
against which the incoming radiation is compared, the pressures of 
the cells can be selected to correspond to the pressure of the 
atmosphere at a particular altitude of interest. Work currently 
underway on the instrument would provide peak responses in the 600 
millibar (mb) and 200 millibar (mb) ranges, which correspond to 
altitudes of 4 km and 12 km. While these altitudes are not partic- 
ularly appropriate to the current need (0-2 km), the possibility 
remains that estimates for lower altitudes could be made by interpre- 
ting data obtained with weighting function peaks at several altitudes. 

CIMATS potentially allows for two other methods of extracting 
vertical information. The first involves the use of the thermal 
(infrared) channel and auxiliary data on the atmospheric temperature 
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and water vapor profiles. This collection of data allows modeling of 
the radiative transfer of the atmosphere, and estimation of the 
vertical profile of pollution. The accuracy of the estimates will be 
a very strong function of the accuracy of the temperature and water 
vapor profile data. The second method applicable to CIMATS would be 
to observe the same pollutant simultaneously in both the solar and 
thermal channels (A and B channels). The solar channel observes 
the total column of the atmosphere from the ground (or cloud tops) to 
the sensor as well as the path from the ground toward the sun. The 
radiation received by the thermal channel will be mainly from the 
altitude to which the instrument is most sensitive (as typified by 
Figure 2-3). Suitable processing of this data could yield information 
on the appropriate altitude range. 

The thermal channels of CIMATS and the MAPS instrument might 
also be able to achieve measures of near-ground pollution by pro- 
cessing the noise in the scene. The signal which reaches the 
instrument is made up of radiation from the atmosphere within the 
f ield-of-view of the sensor (the so-called path radiance) along with 
a contribution from the ground. Generally, the time variation in the 
atmospheric signal is small since transport and diffusion tend to 
make the atmosphere uniform over a relatively large scale. The 
ground signal, however, may vary considerably in time due to the 
nonuniform emissivity of the various components of the surface (land, 
water, vegetation, etc.). So, by preferentially detecting the 
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rapidly varying component of the signal, one can be sure that it has 
been produced by radiation from the ground. This method has been 
developed for a commercially available remote sensor [12]. 

None of the many data processing concepts mentioned in this 
subsection has been successfully used to achieve accurate estimates 
of tropospheric air quality from space. Considerable work would be 
required to establish the most appropriate methods. The use of 
thermal radiation as a source for the sensors requires the simulta- 
neous provision of auxiliary data, particularly water vapor and 
temperature profiles. In addition, information must be available on 
the extent of cloudiness, the distribution of cloudiness within the 
scene and cloud top temperature. Generally, the quality of the data 
produced by the sensor will depend upon the quality of the support 
data. 

Sulfur dioxide and oxidants are particularly important components 
of regional air pollution. Detection of SO^ by both sensors is 
limited by a smaller than acceptable signal to noise ratio. Detec- 
tion of other important pollutant gases (oxidants, hydrocarbons, 
etc.) are not claimed for either instrument. 

A product of the chemistry of sulfur dioxide and nitrogen 
dioxide in their regional transport is a variety of aerosols (both 
sulfates and nitrates). Neither instrument under consideration has 
oeen designed to detect these products of stationary source pollution. 

A last point of concern is the inability of the solar channels 
of CIMATS to perform effectively except under appropriate solar 
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elevation conditions. Figure 2-4 illustrates the yearly variation of 
latitudinal coverage which can be achieved from a 56® orbit if 30® 
solar elevation is required. The Figure shovs limitations in the use 
of CIMATS for various times of the year and also illustrates that 
missions greater than about 20 days are not possible. During that 
period, the motion of the orbit plane relative to the Earth-sun line 
causes considerable changes in the solar elevation along the ground 
track of the orbit. Use of the solar channel of CIMATS in the 
Northern Hemisphere in winter is also precluded above latitudes of 
about 30®. 

In summary, then, the sensors exhibit considerable development 
in their concept and design but appear to be limited in performance, 
particularly from the point of view of obtaining pollution estimates 
in the highly polluted part of the atmosphere near the ground. 
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FIGURE 2-4 

LATITUDES COVERED BY A TYPICAL 56® ORBIT WITH THE REQUIREMENT 
THAT THE SOLAR ELEVATION EXCEEDS 30® 
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3.0 ORBITS 


The characteristics of the orbitCs) chosen for the mission will 
exhibit a considerable influence on data quality. The altitude is 
directly related to the spatial resolution of the observations and 
determines, along with the orbit inclination, the spatial and temporal 
coverage and sampling. Orbit parameters, along with the launch date, 
control the solar elevation along the orbit track v*iich is important 
for the proper operation of CIMATS in its visible light mode. The 
orbit parameters also determine the repeat cycle of the orbit, expressed 
as the number of orbits (or equivalently the number of days) between 
ground tracks which exactly coincide. The repeat cycle controls the 
frequency with which observations are made of a particular site. 

The orbits which might be expected for Shuttle missions will 
range from approximately 250 to A30 kilometers in altitude and from 
28.5 to 57 degrees in inclination. Appropriate selection of values 
from each of these ranges will result in orbits whose repeat cycles 
are an integral numbers of days. (The more general case is that the 
repeat cycle will not be an integral number of days). Examples of 
such orbits are shown in Figures 3-1 through 3-5 and Table “^-l. The 
figures illustrate the ground track which will occur for 1,2, 3, 4, and 
7 day repeat cycles. Each set of ground tracks will repeat after the 
number of days in the orbit's repeat cycle. 

The orbits shown as examples illustrate that as the length of 
the repeat cycle gets longer, the density of ground tracks increases. 

At the same time the frequency of observation decreases. 
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FIGURE 3-2 

GROUND TRACK FOR ORBIT WITH TWO DAY REPEAT CYCLE 
(ahitude » 349 km, inclination « 56^) 







West Longitude 

FIGURE 3-4 

GROUND TRACK FOR ORBIT WITH FOUR-DAY REPEAT CYCLE 
(altitude ■= 277 km, inclination = 57°) 






One of the orbits (Figure 3-1) has an altitude higher than one 
might expect for Shuttle but is included both because it shows the 
control of orbit inclination on latitudinal coverage (no observations 
above 45® occur in the example) and because it has a one-day repeat 
cycle* In effect, this example orbit illustrates one extreme in the 
trade-off between spatial sampling density and frequency of sampling. 
The other extreme, for the proposed seven-day mission duration, is 
shown in Figure 3-5. The orbit has a seven day repeat cycle, about 
four times as much spatial sampling density as the one day repeat 
cycle case but only one seventh the frequency of sampling (since it 
takes seven days to complete the sampling pattern). 

There are two factors which, along with the orbit parameters, 
determine the solar elevation history along the ground track of the 
orbit: the local time of the various observations and the declina- 

tion of the sun as it varies with the seasons. The solar elevation 
is primarily of importance in determining the occasions on which the 
solar wavelength channel of the CIMATS sensor will have adequate 
radiation levels. 

Variation from measurement to measurement of the locel 
time is of significance since it indicates the amount of diurnal 
sampling which can be achieved during the mission. If the mission is 
sufficiently long that each hour of the day and night is sampled many 
times, the diurnal sampling is better than in the case of a sunsyn- 
chronous or nearly sunsynchronous orbit for which all observations at 
a particular point are made at the same local time. 


36 


Seasonal variations in the declination of the sun, and the move** 
ment of the orbit plane with respect to the Earth-Sun line determine 
the range of latitudes for which observations can be made at various 
times of the year. Figure 2-4 illustrates these limits for a 56® 
inclination orbit if 30® solar elevation is to be achieved. The 
Figure shows that missions as long as approximately 20 days can enjoy 
the maximum possible latitude coverage in one hemisphere or the 
other, as long as the launch date is chosen with care. The time of 
day of the launch is also important (in this case it was assumed that 
the first orbit crossed the equator at 12:00 noon). The Figure also 
shows that winter conditions limit the latitudes which can be covered. 
From November to March only latitudes up to 30® can be covered in the 
Northern hemisphere. 

The f ield-of-view of the instrument adds slightly to the area 
covered over and above the sampling pattern of the orbit. For 
example, from an altitude of 400 km, the width of the CIMATS field of 
view (7°), measured in degrees of great circle arc is about 0.4®. 

The spacing of the orbit tracks in the most dense case, a seven-day 
repeat cycle, is about three degrees so the contribution of the f ield- 
of-view is small in comparison. Figure 3-5 illustrates the width of 
the f ield-of-view. 

The orbit analysis model used to develop the figures for this 
section is described in Appendix A. 
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4.0 MODEL FOR STATISTICAL INTERPRETATION 


This section discusses the capability of remote sensors » like 
MAPS and CIMATS, to measure pollution levels when flown on the Space 
Shuttle. This has been done by developing a method of quantifying 
the sensor's performance as a function of the individual instrument's 
design parameters. In particular, the impulse response functions of 
the instruments are derived. These response functions are used along 
with consideration of the fact that the pollution distribution is a 
function of time as well as space. The important design parameters 
considered are: 

1) Velocity wf the ground track of the Space Shuttle, 

2) Resolution ^area on the ground viewed by the instrument), 

3) Temporal response characteristics and response delay, 

4) Time duration of individual observations and, 

5) Spatial and temporal sampling of the region of 
interest. 

These parameters are used, along with the instrument model, to 
develop methods for the determination of the uncertainty associated 
with the observations. 

Characterization of the uncertainty in the observations requires 
information on three factors: 1) the uncertainty in individual 

observations as determined by the performance of the instrument and 
data reduction methods; 2) the uncertainty associated with incomplete 
sampling of the scene in space and time; and 3) uncertainty associated 
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with instrument response to spatial and temporal variability in the 
scene* Identification of instrument uncertainty (both bias and 
random) is relatively straightforward as a result of laboratory 
calibration and field testing. Sampling uncertaint ies are somewhat 
more complex, however, since they involve both the sampling char- 
acteristics provided by orbit, as well as the natural spatial and 
temporal variability of the pollutants* Models of the spatial and 
temporal response of the in«5trument response function were carried 
out to complete the analysis. Calculation of the overall uncertainty 
in the pollution estimates involves a sum of squares calculation 
assuming a linear additive model using the contribut ions from 1), 2) 
and 3) above. Discussion of the overall uncertainties concludes the 
results of this section. 

4.1 Spatial Impulse Response Function 

A number of parameters characterize the capability of a remote 
sensor to respond to the variations in the pollution it is attempting 
to detect. These include: 

velocity of the ground track - V(km/seconds) 

time constant of the instrument* - tiCseconds) 
full angle f ield-of-view - ^2(radians) 

altitude of the spacecraft - A(kilometers) 

footprint of the f ield-of-view - D “ S2A 

Each element which enters the field-of-view acts like a unit 

step function whose amplitude indicates the absolute amount of 


♦Assumes the instrument response to a step function grows as 1-e 
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pollution in that element region. During its time in the field-of- 
view, (D/V), that element produces a response of the form 

After the element has moved outside the f ield-of-view, the response 
approaches zero according to 

-t/r 
e 1 

As the footprint progresses along the orbit path, each successive 
small area interacts with the changed concentration and alters the 
output of the sensor. 

If the pollution level is described as g(x), then the contribu- 
tion from elements in the f ield-of-view at any place x is given 


by [1,2] 


fl (x) * - 


g(x') 1- exp 


(■ •■?)] 


Similarly, the contribution of the elements which have moved 
out of the f ield-of-view is given by 


f2 (x 




dx'. (2) 


The total response of the sensor thus becomes 


,).i j 




exp dx' , (3) 
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Equation 3 is a convolution of two functions and can be rewritten as 
f (x) * ^g(x') hj (x-x') dx' 

where 

(x) is the spatial impulse response function and is given by 


hj (x) 


0 

ij^l-e«p jju(x-D). 


(4) 


U(x-D) here is a unit step function defined as 

x<D 


U(x-D) 


0 

* 


(5) 


x>D 


Equation 4 is simplified further when expressed in terms of 
spatial frequencies. The frequency spectrum (jw^) is the Fourier 
Transform of h^(x); 


90 

H, (ju;,) - ./ h(x) e'H* dx 

i 1 ~00 


1 - e'j-l" 


( 6 ) 


Equation 6 can be plotted as a function of after appropriate 
values of V, D and r have been selected. To obtain a satisfactory 
level of signal against background noise, the time constants of 1 
second and 10 seconds are being currently considered for CIMATS and 
MAPS, repsect ively . The f ields-of-view of the respective instruments 
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are 7® and A. 3® (3). The Shuttle platform velocity, V, is approxi- 
mately 7 km/sec* (2]« 

Equation 6 was evaluated with all the appropriate values of the 
parameters listed above* A plot of with varying spatial 

frequencies is given in Figure 4-1. If the half power bandwidth of 
a system is defined as the frequency where the system response 
function falls to 0*707 of its maximum value, then it is evident 
from the figure that the spatial resolutions of CIMATS and MAPS are 
ISO km and 100 km, respectively* It is interesting to compare these 
figures with the instantaneous f ield-of-view of the two instruments, 
unaffected by the satellite velocity and instrument time constant 
(CIMATS "-50 km and MAPS ^30 km from an altitude of 400 km)* A 
considerable increase in the effective f ield-of-view along the orbit 
track is evident for each instrument (corresponding to poorer resolu- 
tion along the orbit track)* Cross-track resolution is unaffected by 
the dynamic response capabilities of the sensors* 

4*2 Temporal Impulse Response Function 

To determine the temporal impulse response function, the instru- 
ments can be represented by the block diagram as shown in Figure 4-2. 



,(t) 


SAMPLER 


FIGURE 4-2 

BLOCK DIAGRAM REPRESENTATIVE OF A TYPICAL SENSOR 
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FIGURE 4-1 

INSTRUMENT SPATIAL IMPULSE RESPONSE FUNCTION 



The ideal low pass filter in the front end of the sensor limits 
the maximum frequency content in the pollution distribution to be 
observed by the sensor and aids in determining the optimum sampling 
rate. The bandwidth Af of this filter is approximately given by (2] 

Af * 7 — 

where 


is Che electronic time constant of the instrument* 

The output signal g*(t), which is limited to the frequency band 

0 toAf, is sampled at a constant frequency of fo • x , where f >2Af. 

*0 

This sampling rate is required since to transmit a band limited 

signal of duration T. and maximum frequency f , it suffices to send a 

1 m 

finite set of 2f T. independent amplitude samples obtained by sampling 
m 1 

the instantaneous amplitude of the signal at a regular rate of 2f 

m 

samples/second. This is a statement o£ the so-called Nyquist 
criterion. 


g**(t) in Figure 4-2 is a train of impulses given by 

T / 

TT 


5"(t) « ^ 


g (nT ) 6 (t-nT ) 
o o 


(7) 


n»o 


where 6 is the Dirac delta function. 


T is the sampling interval and as discussed above is related 
o 

to Af and the electronic time constant by 

1 




( 8 ) 
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g*it) is a sampled function consisting of short pulses, 

F , of equal duration t at equal intervals Tq. Thus, impulse 
response function of the hold circuit is given by 

ho(t) • U(t) - U(t-T^) for all t^< T . (9) 

w 2 2 


The last filter of the system gives a real time response 
and reflects the time constant the system* Its impulse 

response function is given by 


't/T 


h’ (t) 


1 


( 10 ) 


Finally, the convolution of ho(t) with h'o (t) determines the temporal 
impulse response function. The temporal impulse response function 


h^Ct) is given by 


h^lt) - h(,(t) * ho(t) 

^1-exp 0(t- T^) for all <T . (11) 

The temporal frequency response function is the Fourier Transform 
of h^Ct) and is presented by (jw 2 )* 


00 

^2 “ f h^lt) e 


“ jwjt 


“ JW. T- 
, ■'2 2 

T 

1 jwj (l+jw^Tj) 


dt 


for all Tm <T 


( 12 ) 


The function H (ju;^) was evaluated for • Tq and 0.5 T^, and 
is illustrated in Figure 4-3* The poorer response is evident in 
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those cases where the observation time is smaller than the sampling 
interval . 

4.3 Total System Impulse Response Function 

The pollution distribution is a function of time as well as 
spatial variations. In general, because of meteorology, the pollu- 
tion distribution at any point at time t^, may be influenced by 
the pollution distribution of point at time ^ 

instrument responses in space and time, as derived above, can be 
considered as independent of one another. 

Thus, the total system impulse response function represented 
as h(x,t) is given by 

h(x, t) = h^ (x) • h^Ct) 



Taking the double Fourier Transform of Equation 13, we obtain the 

system bi-frequency impulse response function as 

f f -ju-jX t 

* / f ’.i(x,t) e • e dx dt 

-« -00 

Jwt ^ JlJaT 

1-e ^ ^ ^ 1-e (14) 

jw|DU+jujVTj^) * ^1 * jcujO+jwjTj) . 
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This formulation is applied to instrument performance evaluation in 
the next section. 

4.4 Instrument Performance Evaluation 

Let g^(x,t) be the input pollution distribution at any x, and t, 
and g^(x,t) be the instrument response to this input. To determine 
how accurately the original distribution is reproduced, a criterion 
relating the average power in the input and output pollution distribu- 
tions may be used. The average power of a random process x(t) is 
given by 

E I x^(t)J- 

where 


E is the expectation operator. The difference between the 


average power of g^(x,t) and will be given by 


E-|g^(x,t)| - E ^gjCx.t) j. . 


(15) 


The power spectrum or a spectral density S ( U' , u> ) o r ' h 

gl 1 2 

process g^(x,t) is the Fourier Transform of its autocorrelation: 


00 00 




where 








- 00 -00 


' “ E |gj^(x,t)gj(x+4j, 


(17) 


From Equation 16, using Fourier inversion formala, it follows t^at 


00 OO 




JuJi^ 1 

e ^ e du)j dwj. (18) 


-00 -OO 
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with « o, and * o, the above yields an expression for the aver- 


age power of the input distribution 

00 00 


{*1 ./ j ‘“‘’" 2 * ‘‘“1 


2* 


(19) 


-00 -00 

Then, the difference Z becomes 


® 00 


./ ./ ' ‘,2 / / 


,u)2^ 


-GO “00 


-00 “00 


The spectral density S (u), ,uj.) of the output g^(x,t) is dependent 

§2^2 2 


upon the system impulse response function and the spectral density 
of the input [4]. It is given by 

■ 2 


% ‘“l’“2> 


h( jw > iu)2^ 


Sg^(wi.^2^ 


( 21 ) 


So that 


00 CO 


-03 —CO 
^ X 


1 / / • 2 

~ 7 / H( ju;, >joo^) S (u), ) doj du) 

4^2 / / 12 gj 1 2 12 


( 22 ) 


“ 00 — 

and the percent difference is given by 

/ 00 00 


X r * 100 . 


/ / . . 1 2 
-CO -03 I ju)2) I ‘*'*'1 

1 '" -X -<r 

/ / 


du)« 


2 I (23) 

“X -00 

S (*J^1 >u),^) du) duJ_ 

gj 1 2 12 

\ -X -X 

The difference of Equation 23 can be used as a part of the uncertainty 


analys is « 
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4.5 Use of A Theoretical Input Pollution Distribution Model 


For a demonstration of the instrument error, g^(x,t) can be 


considered a product of two mutually exclusive functions g(x) and 


gj^(t) where g(x) determines the spatial pollution distribution and 


gj(t) give the temporal pollution distribution. 

Furthermore, the spatial distribution of pollutant concentra- 
tions g(x) can be assumed to be a Gaussian distribution. This is 
a convenient assumption, and is used as a basis for air pollution 
models [5]. It is of the form 


g(x) 




exp 




(24) 


where x is the center of the pollution distribution (near the urban 
o 


center), and a is a measure of the plume spread. The range of value 
of a (5 to 30 km) is typical of pollution sources ranging in the 
size from an interchange of several major highways to large cities 
like New York. 


The temporal pollution distribution g^(t) can be modeled on 


the basis of a normal Markoff process with zero mean [6]. 

The autocorrelation function of a Markoff process is given by 

- C U/ 


R(^2) * e 


(25) 


where C is a constant, representing the duration of periods of 
constant conditions. That is, every C seconds the atmosphere is 
assumed to have reached some stable state which persists for C 
seconds. 
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The spectral density of g^(t) is given by the Fourier Transform 


of R(| 2 ) 




Combining Equation 26 with the spectral density of Equation 21, 
we obtain the spectral density of our model of the input: 




f± 


2/0 2 
-X /2a 


-j“>,x 


2 2 

4C*‘+u.‘ 


Finally, substitution of Equations 27 in Equation 23 determines ; 


/ I I mmf • • \ ^ ^ 2C dU)| du>M 

J J \ 

z f. 100 1- 2_® £ (28) 


CD CD 

f f e-"l" 

J J 4C ♦ w 


2 du) j^du) 2 


Equation 28 was evaluated for the MAPS and CIMATS design para- 
metersa The constant C was assumed to take a value of 1000 seconds. 
The results for varying spatial variances in pollution spatial distri- 
butions are shown in Figure 4-4. To determine the effect of system 
temporal response separately, Equation 28 was computed considering 
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! 

only; where is defined in Equation 6. 

These values also are plotted in Figure 4-4. It is evident that in 

the case of both instruments » the spatial response controls the 

total instrument response up to scales approximately the size of the 

instantaneous field-of-view. At much larger scales » up to 300 km, 

the spatial and temporal response contribution is about equal. The | 

instruments operate with ignorably small error at the larger scales. | 

4.6 Uncertainty Due to Sampling Strategies ■; 

One of the key issues in determining the performance of remote 
sensors on-board Shuttle is the determination of the uncertainty 
in assessing pollution distributions because of the sampling distribu- 
tion in space and time. 

The sampling problem arises because the selection of the orbit 
and the inhibiting influence of clouds may make the spatial and 
temporal sampling less than adequate for the purpose of characterizing 
the spatial and temporal variability of the pollution. The adequacy 
of the sampling which is achieved is a function of the variability of 
the pollution in space and time. For a scene exhibiting little 
change with time, only infrequent observations are required. 

For a scene with little spatial variability, a few observations 
at different points may be adequate to obtain statistically meaningful 
averages. The problem at hand, unfortunately, is characterized by a 
high spatial and temporal variability so that sampling limitations 
become important. It is the purpose of this subsection to quantify 
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the sampling by considering both that which can be achieved as well 
as the spatial and temporal variability in the pollution distribution 

The pollutants of interest may exhibit considerable spatial 
and temporal variability* In fact, different pollutants may exhibit 
different patterns of diurnal, seasonal and spatial variability* The 
pollution can therefore be regarded as a stochastic variable, with 
the pollution variations through each day being a record from the 
ensemble of all possible records* Letp(t) represent the ensemble 
mean of P^Ct), the pollution at time of the day, t on the i^^day* 
Then, the average pollution for N days at time t is 

N 

jx(t) • ^ ^ (29) 

i-1 

Note that ^(t) is not independent of time because of the diurnal 
variations, and therefore, the process is a nonstationary stochastic 
process* 

Integrating Equation 29 over the times of the day that instru- 
ments are capable of taking measurements, we obtain the time average 
pollution for N days. This is given by 

T / 

time of the day 

Since Shuttle will be visiting one place only a finite number 
of times in N days. Equation 30 is rewritten as a sum: 


where 


^ .th . 

t . IS the 1 hour 
J 

24N , 

" ‘ zr 


At 




will be determined by the orbit parameters since they control 

the time between passes over the same place. Ideally, one would like 

to keep At as close as possible to the sampling interval given by 

Nyquist criterion, which may be impossible to achieve considering 

large variations in pollution diurnal cycles (see Section 4«2)« For 

example, if the temporal variability of the scene is on the order of 

a day, which is not uncommon, the Nyquist criterion requires that 

samples be made on a scale of fractions of a day. It is unlikely 

that samples can be made with this temporal density. Therefore, 

the estimate of the average pollution at time t ^ only be 

available from a single measurement taken at the hour of the i^^ 

day, P.(t.). Let the variance of this single estimate be a, defined 
I J 


as 


Var 


[Ju.,] 




(32) 


then the estimate of the ensemble average over N days will be given by 


- - E K(t.) 
" j-1 J 


( 33 ) 


and the variance of this pollution estimate, p, will be [7] 
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1 ^ 0 % 

■ ■> I 


i-1, j-1 


where (i-j) is the correlation between the pollution estimates 


and(^(t^) made at different times. 


This analysis can be carried further using the following assump- 


tions: 


e the samples over the period are uniformly distributed over 
the time of the day. 


e the variance in any pollution estimate is v and, 


e the correlation between pollution estimates can be assumed 
to be of exponential form, i.e.. 


(i-j) 


where 


k is the index representing the k repeat cycle, and takes 
values between 0 and n. n is the number of repeat cycles in 
a time period of interest and equals the number of samples 
of a place in that time period. 


o’ is tne correlation coefficient/repeat cycle. If M repre- 
sents the number of days/repeat cycle, and O' is assumed to be the 
correlation coefficient for one day, then 


o' - Mo 


1 


Under these conditions, Equation 34 becomes [8]: 


Var ((i) 



(n-k) e 


-kn' 



(n-l)e* 


-n ♦ e 


-(n-l)«' 


(e"-l)2 


( 36 ) 


Equation 36 can be used to generate a set of curves for different 
values of n, M and o' which illustrate the variance in estimates (in 
units of a) for various mission durations* repeat cycles and pollu- 
tion correlations. Curves developed from this equation are shorn in 
Figures 4-5 through 4-7. The mission durations considered in the 
examples are 7* 15 and 30 days. The data pertaining to 15 and 30 
day missions are included for comparison to illustrate the value of 
missions longer than the expected 7-day Shuttle missions. The cor- 
relation of pollution on two successive days is indicated by p. The 
correlation coefficient* p* for one day is determined from using 
Equation 35* and is given by 

- ln(p) (37) 

Upon examination of Figures 4-5 through 4-7* it is clear that for 
a very high correlation of pollution of two successive days* the 
uncertainty in the estimate of pollution in units ofais not sig- 
nificantly reduced by making a large number of samples. On the 
other hand* for a highly varying pollution scene* the correlation of 
pollution on two successive days is low* and the potential reduction 
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errors due to samrum: in a m oav mission 



in the uncertainty in the estimate of pollution provided by a single 
observation is more pronounced. 

4.7 Application of The Analytical Methods 

The statistical tools developed in the previous sections can 
now be used collectively to determine the data quality %diich might 
result from a real mission. For example, it was found in Section 3 
that a typical Shuttle orbit might have a five-day repeat cycle. 

This information, along with data on the correlation of pollution 
and the probability of cloud cover, can be used to estimate the 
sampling uncertainty. The inability of CIMATS to operate under 
conditions other than relatively high solar elevation angles can 
also be included. 

4.7.1 Pollution Statistics 

An extensive listing of air quality statistics were obtained 
from the Environmental Protection Agency [9]. The data were provided 
in the form of hourly and daily averages of pollution levels for July 
and August 1974. From this data, weekly, bi-weekly and monthly means 
(k) and standard deviations (cr) were computed using conventional sta- 
tistical methods. In addition, the correlations were also computed. 
The data appear in Table 4-1. The correlation computation was 
carried out using the method normally used in regression analysis. 
That is, [10] 
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where all sums range from 1 to n, the number of entries (for example, 
for a 30 day mission, the number of entries will be 30)« Thus, the 
coeff icient p represents the correlation between each measurement and 


the next successive one. 


4.7.2 Use of the Data 


The standard deviation and correlation data of Table 4-1 have 
been used to compute the sampling uncertainty as a function of days 
per repeat cycle. The uncertainty due to sampling is 


Var M 


X 100 


Use of Equation 39 and Table 4-1 result in Table 4-11, which 
is an example of the magnitude of sampling uncertainty which might 
occur in attempting to make estimates of the pollution in the New 
York city area, assuming that an IR sensor, unaffected by solar 
limitations is used. A 30Z chance of cloudiness was assumed, 
based on the results of Appendix B. The sampling uncertainty for 
a mission of the expected duration (7 days) ranges from 20 to 36 
percent, depending on the number of days in the repeat cycle of the 
orbits. However, for a mission of a month, the sampling uncertainty 
is reduced to the range of 10 to 22 percent, depending on the number 
of days per repeat cycle. As described in Section 3, the one-day 
repeat cy^le provides the smallest sampling error. However, the 
orbit which produces a one-day repeat cycle is also the one which 


results in the poorest coverage of the area of interest. 
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TABLE 4-II 


PERCENT SAMPLING UNCERTAINTY IN NEW YORK AS A FUNCTION OF DAYS 
PER REPEAT CYCLE (INCLUDING 50Z CHANCE OF CLOUDINESS) 

PERCENT ERRORS 


MISSION DURATION (days) 


REPEAT CYCLE (days) 

7 

15 

30 

1 

20.14 

14.25 

10.22 

2 

25.16 

17.41 

12.38 

3 

30.32 

20.79 

14.73 

4 

35.77 

23.91 

16.92 

5 

35.77 

26.72 

18.89 

6 

35.77 

29.26 

20.69 

7 

35.77 

31.61 

22.35 


> /vAR (u) 

A 


♦EQUAL TO 


X 100 


Three major elements make up the total uncertainty in estimates 
of pollution. They are: the radiometric uncertainty imposed by the 

instrument performance, uncertainties dictated by the spatial and 
temporal response functions of the instrument, and the sampling 
uncertainty. Under the assumption that these uncertainties are 
independent, and they appear to be, then the data quality is defined 
by a root sum of squares calculation. The total uncertainty is 
given by 

% uncertainty 

where L is the number of observations during the mission 
(T^ : uncertainty in individual estimates 

0 ^ 2 " uncertainty associated with response 

: sampling uncertainty (derived from Table 4-II) , 

I 

Table 4-III illustrates the total system uncertainty for various 

mission lengths and repeat cycles. As in the data of Table 4*11, | 

50 percent cloud cover is assumed. Other assumptions used in the | 

analysis are that each pollution estimate has an uncertainty of 50 

percent and that response uncertainties may be as high as 10 percent 

(based on Figure 4*4). 

An expected result appears immediately; mission duration has 
a significant impact on the uncertainties, the shorter lepeat cycle 
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i , TOTAL UNCERTAINTY (IN PERCENT) IN POLLUTION ESTIMATES 

^ ASSUMING 50Z CHANCE OF CLOUD COVER AND 50Z UNCERTAINTY 

I’ IN INDIVIDUAL POLLUTION ESTIMATES 


MISSION DURATION (days) 


1 

REPEAT CYCLE (days) 

7 

15 


!| 

1 

27.9 

19.4 


i; 

2 

37.1 

25.5 



3 

45.1 

30.9 

4 

i 

4 

52.6 

35.6 

4 

i. 

5 

56.0 

39.8 

4 

? 

6 

59.2 

43.5 


t 

7 

62.3 

47.0 
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orbits exhibiting the best data quality. Conversely, these short 
repeat cycle orbits have the poorest density of sampling of the 
region of interest, as shovn in figures 3“1 through 3~5. The mission 
designer has the choice, then, of attempting to achieve the best pos- 
sible data but for limited locations. Alternatively, he may choose 
to gain the best possible spatial sampling of the region (as illus- 
trated in Figure 3-5) at the expense of data quality at individual 
locations. 

Clearly, a considerable number of assumptions have been included 
in the analysis above (frequency of cloudiness, uncertainty of indi- 
vidual pollution estimates, etc.) It was the object of this particular 
sub-section to illustrate the methods to be used in analyzing sensor/ 
orbit performance. Additional analysis would be required if the 
assumptions were to be modified. 




REFERENCES 


1. L. J. Duncan, E. J. Friedman, E. L. Krltx and E. A. Ward, "An 
Airborne Remote Sensing System for Urban Air Quality," The HITRE 
Corporation, MTR-6601, February 1974. 

2. A. Ghovanlou and E* Friedman, "Evaluation of the Airborne Air 
Pollution Remote Sensor," The MITRE Corporation, MTR-679S, Ma’.ch 
1975. 

3. Information on MAPS derived from a private communication with 
H. Reichle, Langley Research Center, May 27, 1977. Data on 
CIMATS appears in Section 2. 

4. J. B. Thomas, Introduction to Statistical Communication Theory . 
John Wiley and Sons, Inc., pp. 144, 1969. 

5. W. B. Johnson, F. L. Ludwig, W. F. Dabberdt and R. J. Allen 
"An Urban Simulation Model for Carbon Monoxide" Journal of the 
Air Pollution Control Association , Vol 23, No. 6, June 1973. 

6. A. Paopoulis, Probability. Random Variables, and Stochastic 
Processes . McGraw-Hill, pp. 535, 1965. 

7. E. J. Friedman, J. N. Gupta, R. G. Henderson, M. E. Lopez and 
H. J. Sheets, "Evaluation of Earth Radiation Measurements by 
Satellites," The MITRE Corporation, METREK Division, MTR-7410, 
April 1977. 

8. J. S. Bendat and A. G. Piersol, Measurement and Analysis of 
Random Data . John Wiley and Sons, pp. 344, 1958. 

9. Obtained from National Aerometric Data Bank, U.S. Environmental 
Protection Agency, Office of Air Quality Planning and Standards, 
Research Triangle Park, North Carolina. 

10. P. Hoel, Introduction to Mathematical Statistics , pp. 165, 

Wiley, New York, 1962. 


69 


Wa 












5.0 USER INFORMATION NEEDS 

The value of the Space Shuttle as a platform for remote sensing 
of atmospheric pollution will require that the data be of sufficient 
quality to satisfy the needs of potential users. There are basically 
three classes of users for data of this type. They include the 
general public, regulatory agencies and scientific groups. It is 
clear that both the standards for data quality and the format of 
presentation of the data will vary with the users. For example, the 
general public may be satisfied with relatively low quality data 
which depicts the general location and movement of pockets of high 
pollution levels. This data could be further supplemented by 
information of the potential for air stagnation, low level inversions 
or other meteorological phenomena which would tend to produce or 
prolong periods of poor air quality. There would also be, no doubt, 
an interest in determining the transport of pollutants or their 
precursors from city to city. 

Scientific or regulatory users have much more specific require- 
ments for data quality both in terms of accuracy and completeness 
of the sampling. In order to improve our understanding the trans 
port, production and fate of pollutants, detailed studies will be 
required. Tt>e contribution to these studies %rhich can be made by 
remote sensors on spacecraft is limi':ed by their frequency of spatial 
and temporal sampling, frequent interference from clouds, relatively 
poor spatial resolution, limited capacity to resolve layers of the 
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troposphere and temporal gaps between missions. Sections, 2, 3, and 
4 of this report define in detail the source and magnitude of the 
various limitations which would be encountered in attempting to 
apply the instruments of interest, MAPS and CIMATS, to scientific 
or regulatory functions. This section includes a discussion of 
the needs which currently exist for regional air quality data for 
scientific or regulatory purposes and a discussion of the data 
standards which would satisfy the public's demand for general air 
quality information. Section 6 compares these needs with expected 
performance and concludes with a set of recommendations for making 
the suggested missions more able to meet the requirements. 

5.1 Users of Regional Data 

An extensive study of the users and needs for local and regional 
air quality monitoring in the troposphere has recently been completed 
by the National Academy of Sciences [1]. This study provides the 
best and most comprehensive treatment of the topic. Table 5-1 
presents a summary of the results showing which of the various users 
have interests in four basic categories of needs. News media have 
an across-the-board need for information on important features of 
regional air pollution. 

5.2 Measurement Requirements 

The various groups interested in regional air pollution data 
do not have the same lequirements in terms of species measured and 
geographical or temporal scales. The subsections below give 
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a aummary of the compounda recommended for meaaurement or which have 
actually been meaaured in the varioua programa* The data have been 
extracted from a recent NITRE report [2]. 

5.2*1 Air Pollution Speciea to be Monitored 

A aummary of the air pollutanta either recommended for meaaure- 
ment or actually meaaured in varioua regional atudiea ia preaented 
in Table 5~II. The difference in reconunendationa which exiat among 
the referencea derive from their particular goala: 

a National Academy of Sciencea (1] - generaliaeo tropoapheric 
monitoring, 


• Sulfate Regional Experiment [3] - inveatigation of aulfate 
formation from power plant emiaaiona. 

• Stanford Reaearch Inatitute [4] - remote aenaing of regional 
pollutanta, 

• Regional Air Monitoring Sya tern/ Regional Air Pollution Study 
[5] - generalized urban and regional air quality, 

a Los Angeles Reactive Pollutant Program (6] - validation of 
mathematical models of air pollution photochemistry in the 
Los Angeles area. 

a Environmental Protecton Agency/Oxidant Study [7] - inveati- 
gation of rural oxidant levels as related to urban hydro- 
carbons and, 

a Battelle Columbus Leboratoriea [8] - formation and transport 
of ozone in nonurban areas* 

5.2.2 Horizontal and Vertical Meaaurement Spacing 

The horizontal or vertical meaaurement spacing selected for 
regional pollutant measurement must be on a scale which is commen- 
aureate with the scale of geographical changes in pollutant concen- 
trations* Except for precise scientific investigations, the changes 
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occurring in •mall cells or pockets of pollution with diaensions 
messured in meters or tens of iseters ere not generally of concern 
in regional studies* 

Rather, regional studies are directed toward the broad scale 
effect produced over a wide area by an entire urban plume or large 
power plant plume. This fact has lead to the conclusion (4] that, 
"Regionally, horixontal resolution in the range of 10 km to 100 km 
would satisfy most needs*" 

Most regional experiments which have been proposed or conducted 

to date do not achieve the spacing suggested above* The principal 

reason for this is economic* The cost of aircraft or ground-based 

equipment, personnel and operations for such an extensive network is 

usually prohibitive* This sparsity of data has been recognised by the 

designers of the Sulfste Regional Experiment (3), 

"Existing data are not adequate to describe f<illy the chemical 
conversion from SO 2 (sulfur dioxide) to the *ulfates and the 
removal processes of sulfur oxides. There is a great need 
for time-dependent three-dimensional messuresients of many 
aerometric quantities to characterise the evolution and fate of 
the atmospheric SO 2 *" 

It is precisely in this area that spacecraft observations, with 
their ability to provide synoptic data, could be most useful if they 
were able to provide the temporal sampling and spatial resolution 
which is required* 
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5 •2. 3 Temporal Samplir4 ^ 

The time scale specified for regional pollutant monitoring 
should be based on the time scales of the major factors which influ* 
ence regional pollutant concentrations* These include: 

• variations in mobile source emissions caused by the occur* 
rence of "rush hours" in urban areas, 

• variations in power plant or industrial emissions, 

a diurnal and other variations in meteorological variables 
and, 

• chemical and photochemical rates of reaction. 

All four of the above factors have time scales %rhich are less than 

24 hours, indicating a need for observations more frequently than 

daily although a wide range of data*taking rates are useful [4] : 

"Temporal resolution needs in the troposphere vary widely 
and are for the most part poorly defined. For research 
purposes, intervals of from minutes to several tens of minutes 
can be important. At the other extreme., monitoring and control 
standards specify instantaneous peak readings through l*hour 
and 24*hour averages and up to annual averages." 
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6.0 CONCLUSIONS AND RECOMMENDATIONS 


The cepabilitiee of the proposed instruaent/orbit coabinations 
have been aade clear in earlier sections. Individual sites in the 
area of interest will be characterised by estiaates tdiose uncertainty 
will be no smaller than 28 percent if the aissions are of seven dsy 
duration (as currently proposed). Longer aissions will result in 
proportionately smaller uncertainties, as described in Tsble 4-111. 

It is also clear from the orbit traces that the one-day repeat 
cycle, which exhibits the minimum uncertainty, has the poorest cover- 
age of the area of interest. Longer repeat cycles exhibit a denser 
coverage of the area but are characterised by higher uncertainties. 

Data quality is also limited by instrument performance charac- 
teristics. Most important is the fact that no method currently 
exists for directly obtaining pollution data on the loifest altitudes 
of the atmosphere, the region tdiich is most highly polluted. Other 
operating limits, as depicted in Section 2, result in overall perfor- 
mance capability which falls short of the requirements expressed for 
scientific or regulatory missions. Scientific needs require high 
data quality specific to the region of interest with reasonably 
complete spatial sampling. Regulatory uses require high data quality 
over long periods. These performance goals lie outside the range of 
capability of the proposed missions. 

The instrument/orbit combinations do, however, hold the possi- 
bility for providing qualitstive data on the movement of polluted 
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air masaea. Under the conditiona impoaed by the aponaor on inatru- 
ment aelaction and orbit characteriatica, thia qualitative data will 
be aufficient to inform the public of the long-range tranaport of 
polluted air, illuatrate the general relltionahip between meteorology 
and air pollution and, poaaibly, identity major pollution aourcea* 

The reaulta of the atudy alao ahow that miaaion durationa which 

# 

greatly exceed the currently propoaed aeven day flight, will reault 
in data quality %ihich, while not being aufficient tor acientific or 
regulatory purpoaea, may be adequate to provide a reliable picture 
of pollution levela. 

The value of the data to the public will, of courae, be a 
function of ita accuracy and the fidelity with which it deacribea 
the region of intereat* It ia evident, then, that data analyaia 
procedurea which allow determination of pollution levela near the 
ground are greatly needed. The experiment aa deacribed in thia 
report ia arbitrarily limited by the fact that only nadir-viewing 
inatrumenta are conaidered. The incluaion of aome method of croaa- 
track acanning would allow the aelection of abort repeat cycle orbita 
(which minimixe uncertainty in the pollution eatimatea) but atill 
provide good apatial coverage of the region of intereat by acanning 
the f ield-of-view to the areaa which the aatellite doea not fly 
directly over. Thia approach ia not a complete aolution to the 
apatial-temporal aampling problem, however, aince the number of 
aamplea ia atill limited by the duraticn of the miaaion. Scanning 
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off the orbit track will reduce the number of observations of areas 


directly below the spacecraft in order to provide more complete 
spatial sampling between the orbit tracks. The most appropriate 
solution to limited sampling is increasing the duration of the 
mission, as is illustrated in Table 4-111. 
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APPENDIX A 


ORBIT ANALYSIS PROCEDURES 


A.O ORBIT ANALYSIS 


The following sections briefly describe the steps required to 
compute orbit parameters, orbit propagation and solar elevation and 
azimuth angles [1]. A glossary is included along with the logic used 
in computation of the orbit propagation. 

A.l Symbols 

xl * 66058.33128 

x2 = 2tt= 6.283185 

x3 = 398601.2 = G#Me 

x4 = 57.296 = degrees per radians 

x5 = 0.9856473 = angular velocity of the Earth-sun line 

in degrees per day 


G 

Me 

AO 

E 

I 

TO 

MO 

MI 


gravitational constant - 6.67 x 10 j 

Kg sec 

24 

mass of the Earth = 5.983 x 10 Kg 

unperturbed semi-major axis of orbit in km 

eccentricity of orbit 

inclination of orbit in degrees 

unperturbed periods, seconds 

unperturbed mean motion, radians/second 

mean motion, degrees/day, including perturbation due 
to first-order non-sphericity of the Earth 
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;Ci IJLAiiK NOT FILMED 


D1 ■ perigee rate, degreea/day 
D2 ■ node rate, degrees/day 
D3 * perigee step, degrees/mean anomaly step 
D4 » node step, degrees/mean anomaly step 
D5 * an intermediate step in the computations 
D6 * daily synchronous orbit precession 
D7 * yearly synchronous orbit precession 
PO * an intermediate step in the computations 
T1 = anomalistic period 
T2 = nodal period 
Q = repetition factor 
Q1 * synchronous repetition factor 

SO * unperturbed synchronous inclination 

S = synchronous inclination 

D * mean anomaly step 
I * argument of the orbit 

R = longitude between successive ascending nodes 

\ * initial longitude 

o 

^ * solar elevation angle 

4> • latitude 

^ * longitude 

P * solar azimuth 
N ■ orbit number 
b * declination of the sun 
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u) * hour angle of the sun 
n ■ number of days from vernal equinox. 

A. 2 Glossary l2 1 

anomalistic period - the time elapsing between successive 
passages of a satellite through its perigee 

anomaly - the geocentric angle of a point along a satellite 
orbit as measured from perigee 

argument - the geocentric angle of the point measured from 
asrending node in the orbital plane and in the direction of 
motion 

azimuth of the sun - angular distance of the sun from north 
(east is positive) 

declination - angle between the Earth-sun line and the equatorial 
plane 

eccentricity - the ellipticity of an orbit 

hour angle - angular distance of the sun from the meridian 
(or longitude) of the observer 

inclination - the angle between the orbital plane of a 
satellite and the equatorial plane measured at the ascending 
node 

mean anomaly - the geocentric angle which would be swept out 
by a satellite radius vector, since the last passage through 
perigee, if the satellite moved at uniform speed — it is directly 
proportional to the time since passage through perigee 

nodal period - the time elapsing between successive passages 
of a satellite through its ascending node 

node - intersection of the orbit track and the equatorial 
plane of the Earth 

solar elevation - angular elevation of the sun above the 
horizon. 

A. 3 Equat ions 

PO - lAO . U - E^)l^ (1) 
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TO 

MO 

Ml 

D1 

D2 

T1 

T2 


x2 • AO • y AO/x3 
y2/TO 

xl* • (1 - 1.5 sin^ I) 


MO * 1 


[ 


PO 


5 


. 2 


xl « Ml « (2 - 2.5 sin 1) 
PO 

xl » Ml « cos ^ 

PO 

x2 ♦ 180 ♦ 86400 
Ml • TT 

x2 . 180 . 86400 
(Ml Dl) • TT 


( 2 ) 

(3) 

x4 • 86400 (4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


SO * arc cos 


-AO^*^ » PO » x5 
xl • \/x3 . x4 • 86400_ 


.^rc cos 


1 ♦ 

L PO 


cos so 


y 1 - • 


(1 - 1.5 sin^ SO) 


] 


D3 

D4 

D5 

D6 

D7 


xl » (2 - 2.5 sin^ I) » D 
PO 

- xl • cos I • D 


PO 


•[ 


MO • 1 + 


xl 


yr p2 

. V 1 - E • 


(1 - 1.5 sin^ S) 


PO 


xl • D5 • cos S • x4 . 86400 


i»(i 


D6 • 365.2422 


Ml + Dl 


0.7292115061 x 10 • 180 


• 86400 - D2 


Q1 


Ml Dl 
360 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


with Ml and Dl computed at inclination S. (17) 
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A. 4 The Repetition Facto r [3] 


The longitudinal position of the ascending node moves an amount 
each orbit equal to 

R - (q - D2) 3^4 0 0 degrees of latitude (18) 

where 

0 is the rotational velocity of the Earth relative to the stars 
(360.9856473 degrees per day). 

The repetition factor Q, the number of orbits per day, is given 

by: 

Q * 360/R (19) 

Q can be written as 

Q . I ♦ 2 M 360/R (20) 

^ n 

where 

p/n is in its lowejt denominator form and I is an integer. After 
In*»-p orbits, the nodes will have moved 

(In^p)R - n360 (21) 

degrees in longitude and will be back at the starting point (since an 
integer number of circuits of the Earth will have been completed, and 
n is the smallest integer which satisfies Equation 20). The number 
of days required for the orbit to repeat (s repeat cycle) can be 
found by noting that, since Q is the number of orbits per day, 

nQ ■ In-^p (orbits) (22) 

so that the repeat cycle ■ n. 
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That is, n is the number of days required for repeat coverage of a 
particular longitude at the equator (and all other suborbital points 
on that orbit trace)* Thus, any values of p and n represent an orbit 
with an n day repeat cycle* A special case exists if Q is an integer* 
Then, n ■ 1 and a one-day repeat period exists* 

The numoer of orbits within the period of the repeat cycle dete..- 
mines the sequence in which regions are sampled* 

Consider the descending node* If the v»estward displacement of 
the node per oribt is R degrees, then after an integral number of 
orbits, I, the node will have moved IR degrees. If an additional 
fractional orbit occurs each day (that is Q is an integer plus a 
fraction), then definition of the nodal movement per day becomes more 
complex since R can only be measured at the equator. If we define 
the nodal movement per day to be measured to the nearest integral 
orbit, then a Q with fraction < 1/2 will be rounded down to I orbits 
per day. 

Thus, the nodal movement will be 

IR degrees per day for ^ <1/2 . (23) 

n 

If p/n > 1/2 then the orbits per day is rounded up to so 
that the nodal movement is 

(I-^-DR degrees per day for >1/2 . (24) 

Using this method, we see that if p/n < 1/2 for each successive 
day, the movement of the node will move toward the east (since a 
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fraction of an orbit is removed each day). If p/i. > 1/2 the movement 
will be to the west. 

The amount of longitudinal movement of the node is then equal to 

- R for -^ < 1/2 (25) 

n n 

and 

(1 - -^)R for p/n > 1/2 (26) 

n 

noting that the minus sign represents eastward movement. 

A. 5 Orbit Propagation 

The propagation of the orbit in longitude and latitude is con- 
veniently expressed in terms of the argument of the orbit, f. 

o ■ latitude ■ arc sin (sin (180-1) sin H] (27) 

longitude ■ \ ^ arctan [cos (180-1) tan ^1 (28) 

o 

+ £/q ♦ liinSli .3(,o - 180 • P 

where 

n is incremented by 1 each time N-1 exceeds Q 
P - 0 for ^ < 90 

P - 1 for 270 > I >90. 

An additional correction must be made if \ exceeds 180®^: 
if \ <-180® then \ - \ ^ 360® 
if \ >180® then \ - \ - 360® . 

Generally, i is restricted to range from 0 to 360 degrees. As 
a result, the orbit number, N, is increased by one each time r 
reaches 360®. 

*Note: Positive \ is measured westward from the Greenwich Meridian 

and negative X is measured eastward. 
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The orbit propagation is defined by Equations 27 and 28 and the 
initial longitude of the ascending node, , of the orbit. Of course, 
changes from orbit to orbit by 


A . 6 Solar Elevation Angle and Azimuth 

The solar elevation angle, ai is a function of latitude, 4> i 
declination of the sun, b, and hour angle of the sun, u; • 

The declination of the sun refers to its position relative to 
the Earth's equatorial plane and can be expressed as 


6 ^ 23.5 sin 0.9856 n 


where 


n ■ number of days from vernal equinox. 

The hour angle of the sun is defined such that 

^ ■ 0 at noon and equals 90® at 6 p.m. 


The mathematical expression of the solar elevation angle is 


[4,5] 


sin o ■ sin <t> sin 6 cos 6 cos 6 cos 


and the solar azimuth is given by 


sin P 


c o s u s 1 n u 


This expression can be used with the formulae of the previous 
sections to compute the solar elevation angle along the orbit track. 
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The hour angle is given by 


U) * w -f \ 

O O 


\+ l/Q (N-1) (360.9856-D2) (32) 


where 


u)^ is the hour angle at the first equatorial crossing at the 
initial longitude \ . oj is usually taken to be 0 for non- 
synchronous orbits and, or course, equal to the appropriate hour 
angle for synchronous orbits. 

The first term in brackets represents the movement of the sun 

during the time the satellite moves through I degrees of argument 

(0^^060) while the second term represents the movement of the sun 

relative to \ during (N-1) orbits, 
o 

If uj<-i80® thenu)= -*-180® (wis positive west of Greenwich 
Meridian, negative east of Greenwich Meridian). 

To convert to local sun time: 
if w>0t=w/15 

( 180 +u)) 


if w < 0 t 


15 


(positive times are p.m. and negative times are a.m.). 

A. 7 Program Documentation 

Computation of the latitude is straightforward and is described 
by Equation 27. The longitude computation is more complex due to 
the need to keep the range in the conventional limits of 0 to 180 
west of the Greenwich Meridian and 0 to -180 east of the Greenwich 
Meridian, as well as noting the longitudinal movement from orbit to 
orbit and day-to-day. 


i 

i 
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Figure A-1 illustrates the flow and decision process associated 


with the calculation of latitude and longitude. 
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FIGURE A-1 

FLOW DIAGRAM FOR COMPUTATION OF LATITUDE AND LONGITUDE 
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APPENDIX B 


INCIDENCE AND DURATION OF CLOUDINESS IN WASHINGTON, D.C, 

B.O INTRODUCTION 

The basic purpose of this appendix is to analyze the frequency 
of occurrence and duration of cloudy versus clear periods at a 
selected location (Washington, D.C.). This location is assumed to be 
representative of the Eastern U.S, and can, therefore, be used as a 
measure of the frequency with which clouds inhibit operation of 
infrared remote sensors. The average cloudiness statistic, which is 
calculated from measurements at National Weather Service Stations, 
provides no information on these statistics. 

This appendix is presented in two main sections. In the first, 
the clear versus cloudiness statistics which were obtained from the 
National Climatic Center are described. In the second section, the 
ten year annual average data for Washington, D.C. , are analyzed to 
give an example of how the data could be used and what typical 
results indicate. 

B . I Cloud Patterns 

In the temperate climatic zone which is typical of most of the 
United States, clear and cloudy periods can be correlated with the 
passage of high anu low pressure cells and their associated frontal 
systems. This produces cycles of clear and cloudy days, the duration 
of which is somewhat correlated with the speed of movement of the 
fronts and pressure cells. Slower movement produces longer periods 
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of clear and cloudy days and long stagnant periods are usually 
associated with quas i-s tat ionary systems. 

This characteristic of temperate zone weather produces a cycli* 
cal pattern of clear and cloudy days which typically consists of one 
to three cloudy days followed by a similar number of clear days. 
Analysis of data showing the frequency of occurrence and the duration 
of these cloudy and clear periods is necessary as an input to evalua- 
ting the sampling gaps which may occur. 

A complicating factor is introduced when one recognizes that 
partly cloudy conditions do occur and, in fact, are probably the most 
usual case over much of the United States. Thus, the analysis would 
consider a selected number of categories of partly cloudy conditions 
in addition to completely clear and completely cloudy cases. The 
other factor which could be considered is that clear, cloudy or 
partly cloudy conditions do not occur in convenient one day incre- 
ments but may change from hour to hour or even more frequently in 
some cases. At this level of complexity, analysis of the actual hour 
by hour reports of cloud amounts at the location of interest would be 
necessary. 

In order to achieve the basic purpose of this appendix, i.e., 
demonstrating the statistics of the cloudiness data, it was decided 
to obtain the necessary data to perform an analysis as discussed 
below. The basic statistic of interest would be the frequency of 
occurrence of cloudy periods of n days duration which are preceded by 
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a clear period of m days duration. The National Climatic Center was 


contacted to provide such data. Since the data are nonstandard, it 
was necessary to prepare a specific computer program which would 
generate the data from standard observations. This program had the 
following specifications: 

1. Based on standard National Weather Service hourly (during day- 
light) observations of total cloud cover, each day at each 
station of interest was categorized as either clear or cloudy by 
determining the arithmetic mean of the total cloud cover for all 
hours from sunrise to sunset. Days with five-tenths or less 
mean total cloud cover were categorized as clear and days with 
six-tenths or more total cloud cover were categorized as cloudy. 
The simple assumption is therefore made that the amount of 
sunshine available during periods when total cloud cover is 
between six-tenths and nine-tenths is offset by the amount of 
sunshine which is lost during periods when total cloud is 
between one-tenth and five-tenths. 

2. Cycles of clear and cloudy days were determined and counts were 
made of the number of occurrences n cloudy days preceded by m 
clear days, n and m were presented in eleven categories each (1 
to 10 and '^lO). 

3. Observations for the ten-year period from 1965 through 1974 were 
analyzed . 
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The results were presented by individual season and annually for 

each of the ten years and for the total of the ten-year period. 

Seasons were defined as: 

Winter - December to February 
Spring - March to May 
Summer - June to August 
Fall - September to November 

Those episodes which extended over two seasons or two years were 
arbitrarily categorised into the time slot corresponding to the 
first clear day of each clear/cloudy cycle. 

A sample of this data is shown in Table B-I which presents the 
average values of cloudiness for Washington, D.C. for the period 1965 
to 1974. This particular data set was used in the analysis presented 
in Sect ion 4. 

The table illustrates that the most probable pattern of cloudi- 
ness is one in which a clear day is followed by a cloudy day. Also 
highly likely is the case that a clear day is followed by two cloudy 
days. Extended periods of cloudy or clear conditions become less 
likely as the time interval is extended. It is also interesting to 
note that the total number of cloudy days in the period of study was 
2157 while the number of clear days was 1462. 

Inspection ot the totals for cloudy periods illustrates that of 
the 721 periods which occurred (of varying durations) over 70 percent 
were of three davs or less. 
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